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Summary 

A  random-coefficient  linear  function  of  two  independent  exponential 
variables  yielding  a  third  exponential  variable  is  used  in  the  construction 
of  simple,  dependent  pairs  of  exponential  variables.  By  employing  antithetic 
exponential  variables,  the  constructions  are  developed  to  encompass  negative 
dependency.  By  employing  negative  exponentiation,  the  constructions  yield 
simple  multiplicative-based  models  for  dependent  uniform  pairs.  The  ranges 
of  dependency  allowable  in  the  models  are  assessed  by  correlation  calcula¬ 
tions,  both  of  the  product  moment  and  Spearman  types;  broad  ranges  within 
the  theoretically  allowable  ranges  are  found.  Because  of  their  simplicity, 
all  models  are  particularly  suitable  for  simulation  and  are  free  of  point 
and  line  concentrations  of  values. 
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1 .  BACKGROUND 


This  paper  introduces  some  models  for  pairs  of  dependent  exponential 
randan  variables  constructed  as  random-coefficient  linear  functions  of  pairs 
of  independent  exponential  random  variables.  The  transformations  used 
occurred  first  in  connection  with  time  series  models  in  exponential  variables 
(see  Lawrance  and  Lewis  [1981]  for  details  of  this  work).  These  models 
indicate  that  a  specific  random-coefficient  linear  function  of  independent 
exponential  variables  leads  to  further  univariate  and  multivariate  exponen¬ 
tial  variables,  somewhat  analogously  to  linear  functions  of  independent 
Gaussian  variables  leading  to  further  univariate  and  multivariate  Gaussian 
variables.  Negative  exponentiation  of  the  transformation  then  leads  directly 
to  a  random-power  multiplicative  transformation,  and  it  is  found  that  this 
can  similarly  be  employed  ir.  constructing  pairs  of  dependent  uniform  random 
variables  from  pairs  of  independent  uniform  variables.  Both  transformations 
are  particularly  suited  to  simulation  work,  the  pairwise  uniform  giving  a 
base  from  which  to  transform  to  many  other  bivariate  distributions. 

This  paper  then  presents  results  on  both  dependent  exponential  pairs 
and  dependent  uniform  pairs.  The  transformations  used  are  necessarily 
limited  to  situations  of  positive  dependency,  but  broadening  to  negative 
dependency  is  achieved  by  employing  simple  antithetic  operations  in  the 
constructions.  The  distributions  given  are  free  from  point  and  linear  con¬ 
centrations  of  probability.  If  these  are  required,  as  is  sometimes  the  case 
in  reliability  studies,  the  concentrations  of  probability  can  be  put  in  by 
using  mixtures. 

There  are,  of  course,  many  earlier  bivariate  exponential  and  uniform 
distributions;  for  instance,  the  work  of  Marshall  and  01 k in  [1967]  which 
involves  choosing  the -minimum  of  two  exponentials,  and  which  can  easily  be 
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simulated,  but  which  always  has  a  line-concentration  of  probability.  Downton 
[1970]  and  Hawkes  [1972]  developed  bivariate  exponentials  employing  the  idea 
of  bivariate  geometric  compounding  of  independent  exponentials;  in  simulation 
this  would  require  simulation  of  bivariate  geometries.  Gaver  [1972]  proposed 
a  bivariate  exponential  specifically  for  negative  dependency,  also  allowed 
in  Hawkes'  extension,  and  which  in  simulation  would  need  both  the  operations 
of  choosing  minima  and  ordinary  geometric  compounding;  further,  it  is 
intrinsically  non-reversible,  unlike  most  bivariate  exponential  models. 

Other  types  of  bivariate  exponentials  could  clearly  be  obtained  via  the 
inverse  probabil ity- integral  transform  of  other  standard  bivariate  distribu¬ 
tions;  bivariate  uniforms  in  this  category  are  specifically  considered  by 
Barnett  [1980].  All  these  methods  involve  a  fair  amount  of  complication, 
particularly  the  inverse  probability  integral  transform,  and  limited  flexi¬ 
bility.  The  methods  proposed  in  this  paper  seem  simpler  than  all  previous 
methods  at  least  in  the  major  respect  that  each  bivariate  dependent  pair  is 
obtained  by  simple  random-coefficient  linear  transformation  of  a  correspond¬ 
ing  independent  pair;  furthermore  they  are  flexible  in  their  dependency 
properties,  and  easy  to  simulate. 

In  summary,  unless  very  specific  modelling  or  structural  detail  is 
required.e.g.  the  linear  regression  of  Downton's  model,  the  present  models 
offer  the  following  advantages: 

(i)  The  constructions  are  very  simple. 

(ii)  They  require  only  two  independent  exponential  random  variables 
and  one  or  two  binary  random  variables  for  their  construction. 

( 1 i 1 )  There  is  a  broad  range  of  the  attainable  dependency,  as  measured 
linearly  by  the  product-moment  correlation  or  monotonically  by 
the  Spearman  correlation. 


(iv)  The  models  are  analytically  tractable  so  that,  in  some  cases, 
closed  form  joint  probability  density  functions  and  regressions 
can  be  obtained. 

The  scheme  of  the  following  sections  of  the  paper  is  to  examine  the 
properties  of  the  very  simplest,  one  parameter,  models  (EP1  and  UP1)  in 
detail.  The  properties  include  product-moment  and  Spearman-type  correla¬ 
tions,  joint  probability  functions  and  regressions.  These  models  have  a 
broad  but,  positive  correlation  range  which  is  then  extended  to  negative 
correlation  by  using  antithetic  pairs.  In  later  sections  models  with  more 
parameters  which  can  attain  correlations  right  up  to  the  theoretical  bounds 
are  introduced,  but  developed  in  less  detail.  The  most  general  of  these, 
the  EP+  and  EP-  models,  allow  the  complete  range  of  possible  correlations 
for  bivariate  exponential  pairs  and  the  possibility  of  acconmodating  other 
structural  details;  they  include  all  the  other  exponential  models.  Parallel 
development  is  given  for  bivariate  uniform  pairs. 
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2.  ONE  PARAMETER  POSITIVELY  DEPENDENT  EXPONENTIAL  MODEL:  EP1  MODEL 
The  random  coefficient  linear  function  referred  to  in  Section  1  is 
central  to  this  work  and  is  presented  first.  Let  a  pair  of  independent, 
identically  and  exponentially  distributed  random  variables  be  denoted  by 
(Ep  E2 ) ;  for  simplicity,  unit  means  are  assumed  throughout.  Let  I  be 
a  binary  random  variable,  independent  of  (E^,  Eg),  with  distribution  given 
by  P(I  =  0)  =  e  ,  P(I  =  1 )  =  1  -  6  .  Then  the  random  variable  X  ,  given 

by 

X  =  BE1  +  IE2,  (2.1) 

has  an  exponential  distribution  with  unit  mean;  this  can  be  seen  by  taking 
moment  generating  transforms, 

<f>x(t)  =  E{exp(-tX)}  =  1  (0  +  (1-0)  y-^}  =  yyjr  (2.2) 

The  result  can  be  extended  by  taking  Eg  to  be  a  random-coefficient  linear 
function  of  a  furthe*  pair  of  exponential  random  variables  independent  of 
E^  .  However,  (2.1)  is  all  that  is  required  here  to  construct  bivariate 
pairs  of  exponential  random  variables;  the  result  (2.1)  is  implicit  in  the 
construction  of  exponential  time  series  models  evolving  out  of  the  EAR(l) 
model  of  Gaver  and  Lewis  [1980]. 

The  EP1  Model 

The  first  and  simplest  bivariate  exponential  distribution  of  the  paper 
is  obtained  by  interchanging  E^  and  Eg  in  (2.1)  to  give  the  second 
member  of  a  pair,  and  hence  produce 
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X1  =  6E1  +  E2  , 


(2.3) 


(EP1) 

%2  ~  ^2  +  " 

These  will  be  called  the  EP1(0)  pair;  the  pair  is  clearly  dependent  because 
of  the  comnon  E^,  E2  and  I  .  There  are  clear  general izations  obtained  by 
using  different  B's  and  I's  for  X-j  and  X2  which  will  be  considered  later. 

The  dependence  of  (X^  X2)  can  be  assessed  by  their  (product  moment) 
correlation  coefficient  as 

Corr(Xr  X2)  =  E{(6E1  +  IE2 ) (BE2  +  IE]  )}  -  E(gE]  +  IE2)E(0E2  +  ) 

=  30(1-0)  ,  (2.4) 

which  is  sketched  in  Figure  1.  The  correlation  thus  has  the  non-negative 
range  (0,  0.75)  with  maximum  value  0.75  at  0  =  0.5  and  symmetry  about 
0  =  0.5.  The  product  moment  correlation  emphasizes  linear  dependency, 
taking  maximal  values  of  +  1  if  and  only  if  the  variables  are  linearly 
related  which  is  not  the  case  here;  a  more  general  and  more  appropriate 
measure  of  non-linear  dependency,  Spearman's  correlation,  will  be  considered 
in  Section  3. 
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Figure  1.  Correlation  functions  for  the  EP1(6)  model 
and  the  UPl(g)  model  as  functions  of  £;  both  models 
exhibit  only  positive  dependence. 

The  joint  probability  density  of  X-|  and  X2  can  be  obtained  by 
considering  the  probability  element  for  X-j  and  X2  in  terms  of  8E-|  +  IE2 

and  BE2  +  IE ^  .  On  letting  I  take  its  values  1  and  0  with  probabilities 
1  -  3  and  6  this  becomes 

(1-BjPUBE,  +  E?)  €(x,,  x,  +  dx,},  (BE2  +  E1 )  €(x2>  x2  +  dx2)} 

1  ^  (2.5 

+  8P(BE1  €(xr  x1  +  dx^PigEg  €(x2,  x2  +  dx2)}  . 

The  first  line  of  (2.5)  can  be  expressed  as  a  probability  element  statement 
in  terms  of  E-|  and  E2  separately,  and  using  the  Jacobian,  (1-32),  of 
the  transformation  from  E-j  and  E2  leads  to  the  joint  pdf  of 
X^  and  X2  as 
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(2.6) 


fXrX2(xrx2)  =  (l-6)l(x1  >  gx2,  x2  >  3x1)f(1_B2)E  (x2-Bxi) 

x  f(l-62)E2(xreX2)(1'02)  +  6f3E^Xl^f6E2^x2^  » 

where  1(  )  is  an  indicator  function  taking  the  value  1  if  x-,  >  0x2  and 
x2  >  6x-j  and  0  otherwise,  and  f  denotes  the  density  of  the  exponential 
variable  in  its  suffix.  Simplifying  (2.6)  then  gives  the  desired  result 

fX  X  (xi’x2*  =  1(8x2  <  X1  <  exp{-(x1+x2)/(1+&) > 

1>2  (2.7) 

+  0”^  exp{-(x^+x2)/B>  ,  0  <_  0  £  1 ;  x^  ,x2  £  0  . 

The  joint  density  thus  has  a  wedge  shaped  sector  of  relatively  high  density 
(px2  <  x-j  <  B ”^x2)  superimposed  over  an  independent  bivariate  exponential 
density.  This  function  is  isometrically  plotted  in  Figure  2  for  B  =  0.5. 

For  low  values  of  B  the  raised  section  spreads  out  over  the  quadrant  and 


Figure  2.  Isometric  plot  of  the  joint  probability  density 
function  of  the  positively  dependent  EP(b)  model  for  B  =  0.5. 

This  value  of  B  gives  the  maximum  attainable  correlation  of  0.75; 
the  figure  has  been  truncated  for  both  variables  at  0.5. 
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leads  to  independence  at  8  =  0;  for  high  values  of  8  the  raised  section 
diminishes  towards  nothing  and  again  gives  independence  at  8=1. 

This  joint  density  function  could  be  used  to  obtain  a  maximum  likelihood 
estimate  of  8  from  a  bivariate  sample;  note  that  a  unique  moment  estimator 
can  not  be  obtained  from  the  correlation  (2.2). 


The  joint  moment  generating  function  of  the  (X1,X2)  pair,  <j>x  x  (t,,t2),  i 
easily  obtained  via  the  EP1  defining  equation  (2.3);  thus 


$X  x  (^i ,t2)  =  E (exp(-ti Xi  ~  t2X2) } 

=  E(exp[-(8t-|  +  It2iE1  -  (8t2  +  It1  )E23> 

=  8(i  +  et1  )-1  (1  +  8t2)_1 

+  (1-6)(1  +  Bt1  +  t2)_1(l  v  t1  +  8t2)_1  .  (2.8) 


Differentiating  with  respect  to  t2  at  t2  =  0  and  inverting  with  respect  to 
t1  allows  the  regression  of  X2  on  X-|  to  be  obtained  as 

E(X2|X1  =  x)  =  ex  +  8  +  1  -  (1  +  x/8)exp{-(8_1  -  1 )x>  .  (2.9) 


This  conditional  expectation  will  be  fairly  linear  except  for  small  x  ,  and 
the  asymptotic  conditional  variance  will  be  linear  for  large  x  .  Note  from 
the  joint  moment  generating  function  (2.8)  that  the  pair  (X^,X2)  is  reversible. 


‘J’X^ (VV  =  *x1  .X^VV 


in  particular,  E(X^|X2  =  x)  =  E(X2|X-|  =  x)  . 


An  aspect  of  this  reversibility  is  that  P{X2  >  X^ }  =  0.5  which  can  be  shown 


directly  from  the  definition  (2.3). 
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3.  THE  NEGATIVE  EXPONENTIAL  TRANSFORMATION  TO  UNIFORMITY:  UP1  MODEL 

It  is  well  known  that  for  an  exponential  random  variable  X  of  unit 
-X  -X 

mean,  both  e  and  1-e  are  uniform  random  variables,  distributed  over 
(0,1).  These  results  are  of  double  interest  here.  Firstly,  application  of 
the  transform  to  both  components  of  a  dependent  exponential  pair  will  lead 
to  a  pair  of  dependent  uniforms,  and  secondly,  which  may  be  more  important 
to  stress,  the  product-moment  correlation  in  this  pair  of  uniforms  is  a 
population  analogue  of  Spearman's  rank  correlation  for  the  original  dependent 
exponential  pair.  This  was  demonstrated  generally  by  Kruskal  £1958];  he 
showed  that  the  product-moment  correlation  in  the  transformed  pair  of 
dependent  uniforms,  could  be  linearly  related  to  the  concordance  (quadrant) 
probability 

PTX-,  -  X2)(Y]  -  Y3)  >  0)  , 

where  (X .  ,Y.. :  i  =  1,2,3)  are  three  independent  pairs  with  the  untransformed 
bivariate  distribution.  Estimation  of  this  probability  was  based  on  con¬ 
structing  the  empirical  joint  distribution  of  an  observed  sample,  and  led 
Kruskal  to  estimating  the  required  linear  function  of  the  concordance  prob¬ 
ability  by  Spearman's  correlation.  Thus,  we  have  the  useful  idea  that 
product-moment  correlation  in  a  bivariate  distribution  after  its  transforma¬ 
tion  to  uniform  marginals  represents  a  population  analogue  of  Spearman's 
correlation  for  the  untransformed  bivariate  distribution.  As  such,  it  is  a 
suitable  measure  of  monotonic  dependency  for  a  bivariate  distribution.  In 
particular,  the  Spearman  correlation  is  a  suitable  measure  for  dependency  in 
bivariate  non-Gaussian  distributions  because  of  their  intrinsic  non-linearity, 
and  the  consequent  limitations  of  the  product  moment  correlation. 
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General  Results  for  the  Transformation 

In  the  case  of  an  exponential  pair  (Xj,X2)  of  unit  means,  the  preferred 
transformation  to  a  uniform  pair  is  given  by 

^  =  1  -  exp(-X1),  Y2  =  1  -  exp(-X2)  .  (3.1) 


The  correlation  and  regression  for  this  pair  are  easily  obtained  from  the 
joint  moment  generating  function,  $x  x  (t|,t2)  =  E{exp(-t-jX1  -  t^)} 
of  (X-j,X2)  .  First  for  the  correlation  we  have 

E(Y1Y2>  =  E{(l-exp(-X1)(l-exp(-X2)}  =  *x  x  (1,1)  ,  (3.2) 

and  thus  the  correlation  result 


Corr(Y1,Y2)  =  12*x  x  (M)  -  3  • 


(3.3) 


For  the  regression  of  Y^  on  Y2  ,  first  write  it  in  terms  of  X1  and  X2  as 

E(Y1|Y2  »  y)  =  E(l~exp(-X-j)|X2  =  -  log(l-y)}  .  (3.4) 

Now  setting  x  =  -  log(l-y)  and  denoting  the  conditional  p.d.f.  of  X1 ( X2 
by  jx  (x,|x)  ,  the  regression  is 

EtY1|Y2  -  l-exp(-x)]  =  /  [l-exp(-x^ )]fx  jX  (x^|x)dx1  .  (3.5) 

x-|=0  12 

The  required  result  is  obtained  after  multiplying  both  sides  by  exp(-tx) 
and  integrating  x  over  (Q,~),  so  leading  to 
« 

f  exp(-tx)E[Y.  |  Y  =  1  -  exp(-x)]dx  *  t  -  <{>¥  y  (l,t-l)  .  (3.6) 

o  1  L  Al’*2 
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i 

Thus,  the  regression  of  the  uniforms  follows  from  inverting  the  joint  moment 
generating  4>v  Y  (l,t-l)  of  the  exponential  pair  with  respect  to  t  as  a 

A-J  ,  A2 

function  of  x  ,  and  then  replacing  x  by  -  log(l-y)  . 

The  UP!  Model 

The  negative  exponential  transformation  to  uniformity  is  now  exemplified 
using  the  EP1(6)  model.  To  obtain  a  pair  of  independent  uniform  random 
variables  from  a  pair  of  independent  exponentials  it  is  simplest  here  to  use 
the  monotonic  decreasing  version  of  (3.1),  that  is 

U-,  ~  exp(-E1 )  ,  U2  =  exp(-E2)  .  (3.7) 

It  follows  then  that  the  dependent  uniform  pair  derived  from  the  EPl(e) 
exponential  model  via  (3.1)  is  given  by 

(UP1)  Y1  =  1-uf  u\  ,  Y2  =  l-u{  U8  ,  (3.8) 

where  I  retains  its  definition  from  Section  2.  The  random  power  multipli¬ 
cative  aspect  of  the  model  is  evident.  The  pair  (Y-|,Y2)  of  (3.8)  will  be 
termed  the  UP1(8)  model;  the  correlation  of  this  uniform  pair  follows  imme¬ 
diately  from  the  joint  moment  generating  function  of  (X-j  ,X2 )  as  given  at 
(2.8),  or  by  direct  calculation,  and  is 

Corr(Yi,Yj  -  38(1-8)  - 8*-7— - —  .  (3.9) 

1  2  ( 1 +3) 2 ( 2+6) 2 

This  function  of  8  is  illustrated  in  Figure  1;  its  maximum  value  is  0.72 
for  8  *  0.32.  The  joint  density  of  Y]  and  Y2  may  be  obtained  by  considering 
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(3.10) 


P(uf  Ug  £.  y-j »  Ug  U-j  £  y2) 

■  6P(U®  <  y1)P(U^  <  y2)  +  (l-B)P(uf  Ug  <  »  U]  Ug  <  y2)  . 


Density  for  the  last  probability  is  restricted  to  the  leaf  shaped  region 
(y^e  <  <  yg,  0  <y 2 <  1).  and  the  final  result  takes  the  form 

fYrY2(yiy2*  3  1(y2/eiy1  <  ygKi+er^y^r6^1'^  +  B‘1(y1y2)^1"^/B 

(3.11) 

(0  i-y^yg  <1;  0  <  8  <  1)  . 

The  density  is  isometrically  plotted  in  Figure  3  for  b  =  0.32,  the  value 
giving  maximum  product  moment  correlation. 


Figure  3.  Isometric  plot  of  the  joint  probability  density 
function  of  the  positively  dependent  UP1(B)  model  for  3  =  0.32  . 
This  value  of  8  gives  the  maximum  attainable  correlation  of  0.75. 


* 
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Regression  for  the  UP!  model  is  most  easily  obtained  via  (3.6)  using 
the  joint  moment  generating  function  of  and  X2  given  at  (2.8)  and  takes 
the  form 

E(Y2| Y1  =  y) 

=  1  -  (l+6)-1(l-y)(1"6)/3  +  (2+6)'1(l-y)(2‘e)/B  -  (2+8)_1(l-y)B  .  (3.12) 

The  regression  has  initial  value  0/(1 -B)  which  is  always  less  than  0.5  and 
curves  upwards  towards  1.0  in  a  near  monotonic  concave  way.  For  6  =  0  and 
6=1  it  has  the  value  1/2  which  is  correct  since  Y2  and  Y-j  are  then 
independent  and  E(Y2)  *  1/2.  Note  that  the  UP1  pair  is,  like  the  EP1  pair, 
reversible,  in  particular  E(Y2|Y^  =  y)  =  E(Y||Y2  B  y)  .  Also 
PtU-j  >  U2>  =  0.5;  this  can  be  seen  directly  from  the  definition  (3.8). 


13 


4.  NEGATIVE  DEPENDENCY  USING  AN  ANTITHETIC  APPROACH:  EP2,  UP2  MODELS 
The  antithetic  transformation  of  any  random  variable  X  is  found  by 
first  expressing  it  as  a  function  (the  Inverse  probability  integral)  of  a 
uniform  (0,1)  variable  U  ,  and  then  by  taking  the  same  function  of  1-U.  For 
any  nonnegative  random  variable  X  ,  the  antithetically  transformed  variable 
a(X)  has  the  same  marginal  distribution  as  X  and  is  maximally  negatively 
correlated  with  X  ;  note  that  the  Spearman  correlation  of  X  and  a(X)  is 
minus  one,  but  their  product  moment  correlation  will  not  necessarily  be  as 
negative  as  this.  In  the  exponential  case, 

a(X)  =  -  log(l-exp(-X))  (4.1) 

and  the  maximum  negative  product  moment  correlation  is  -  0.6449,  as  is  well 
known  after  Moran  [1967].  As  a  pair  of  exponential  random  variables,  X  and 
a(X)  lie  on  the  curve 


e-x  +  e-a(x)  ,  , 

and  so  are  highly  degenerate  in  the  sense  that  X  completely  determines  a(X)  . 

The  EP2  Model 

The  aim  here  is  to  construct  negatively  dependent  exponential  pairs 
which  are  not,  like  the  antithetic  pair,  partially  or  completely  degenerate; 
the  random-coefficient  linear  transformation  (2.1)  allows  us  to  do  this  while 
still  making  use  of  antithetic  pairs.  The  general  approach  is  illustrated 
with  respect  to  the  EP1  model;  in  (2.3)  the  E^  and  E2  \  are  replaced  by 
their  antithetic  versions  a(Ej),  a(E2),  and  I  In  X-|  and  X2  Is  replaced  by 
1^  and  I2  respectively,  these  variables  having  the  same  marginal  distribution 
as  I  .  Thus  the  suggested  negatively  dependent  model  has  the  form  ? 
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X.j  =  0Ei  +  1 1  ^2  * 

( EP2 )  (4.2) 

and  the  expression  for  its  correlation  is  given  by 

Corr(XrX2)  =  20(1-B)(-O.6449)  +  Cov(IrI2)  .  (4.3) 

This  suggests  that  maximum  negative  correlation  of  X^  and  X2  requires 
maximum  negative  covariance  of  I-j  and  I2  .  For  fixed  3  ,  this  will  be 
achieved  for  minimum  E(IpI2),  or  equivalently  minimum  p  =  P(I-j  =1,  I2  =  1)  . 
Consideration  of  the  general  joint  distribution  of  I,,  I2,  that  is 


!1  * 

1 

0 

I2. 1 

p 

1  -0— p 

1-B 

0 

i-e-p 

2B-l+p 

8 

i-b 

B 

1 

(4.4) 


shows  that,  since  the  probabilities  must  be  non-negative,  the  (0,0)  term 
constrains  p  to  the  minimum  value  of  l-2s  if  8  <_  h  and  0  if  0  >_  %  .  The 
implied  joint  distribution  is,  not  surprisingly,  the  antithetic  one,  and  it 
has  covariance  -  B2  for  0  <  6  <  H  and  -(1-6) 2  for  H  <  6  1  1  .  Thus,  going 
back  to  (4.2),  the  model  *••' th  1^  =  I,  I2  =  a(I)  is  of  most  interest  in 
connection  with  negative  dependence  and  will  be  called  the  EP2(6)  model;  its 
correlation  is  given  by 


Corr(X.|  ,X2) 


'  -2B(1-B)(0.6449)  -  62 
-200-6X0.6449)  -  ( T-0) 2 


0  <  6  <  h  , 
h  <  3  <  1  . 


(4.5) 
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This  function  is  sketched  in  Figure  4.  The  minimum  value  is  -0.5724  at 
3  =  0.5  which,  considering  the  model  is  never  degenerate,  is  fairly  satis¬ 
factory.  (The  Gaver  model  [Gaver,  1972]  is  also  not  degenerate  and  has  a 


Figure  4.  Correlation  functions  for  the  EP2(g)  model 
and  the  UP2(3)  model  as  functions  of  6;  both  models 
exhibit  negative  dependence. 

minimum  correlation  of  -0.5).  A  simulation  of  EP2(3)  pairs  at  the  B  =  0.5 
value  is  given  in  Figure  5;  the  scatter  is  best  described  as  a  non-degenerate 
emulation  of  the  antithetic  concentration  along  the  curve  e~x  +  e”y  =  1  . 
Weakening  of  the  dependency  as  the  parameter  3  approaches  0  or  1  fades 
this  characteristic  scatter  into  the  blander  independent  exponential  picture. 
The  scatter  in  Figure  5  should  be  compared  to  the  isometric  plot  of  the  joint 
density  in  Figure  3. 

The  joint  moment  generating  function  of  X-j  and  X2  requires  the  joint 
moment  generating  function  of  Ei  and  a(E.)  which  is  straightforwardly 
evaluated  as 
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Figure  5.  Scatter  plot  of  50  simulated  pairs  of  observations 
from  the  negatively  dependent  EP2(B)  model  for  B  =  0.5:  the  solid 
line  represents  the  curve  relating  the  antithetic  bivariate 
exponential  pairs. 


1  t,  t9 

<p(t,,t2)  =  /  x  1  ( 1  -x )  ^dx  ,  (4.6) 

c  o 

which  is  the  Beta  function  B(1  +  tj ,  1  +  t2)  .  Also  required  is  the  joint 
distribution  p^  =  P{I  =  i,  a ( I )  =  j),  i,  j  =  1,0,  available  as  remarked 
previously  following  (4.4);  hence  the  result,  via  (4.2), 

■  P114>(3t1  ,t2)<p(t1 ,3t2)  +  p^UV^Ml+B^)"1  (4.7) 

+  P014>(t1,6t2)(l+Bt1)_1  +  pQQd+Bt^’^l+Btg)"1  . 
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Actually  the  model  (4.2)  works  for  any  independent  bivariate  exponential 

I  I 

pairs  ( E-j  ,E^ )  and  (£2^2)  •  However,  the  antithetic  pair  is  simplest  to 
generate  and  gives  the  greatest  negative  correlation  in  the  (XpX2)  pair. 

The  UP2  Model 

Thus  far  we  have  the  EP2  model  giving  a  broad  range  of  negative  depen¬ 
dency  for  exponential  pairs;  this  is  now  transformed  to  the  corresponding 
uniform  pair  giving  the  UP2  model.  By  negative  exponentiation  of  (4.2)  and 
use  of  (3.1)  and  (3.6)  we  have 

(“K)  V,  -  1  -  V2  '  1  -  0-U,)a(I)(l-U2)B  .  (4.8) 


The  correlation  of  and  Y2  may  be  evaluated  directly  from  (4.8)  or  by  use 
of  (4.7)  and  (3.3)  and  gives  the  result 


CorrfYpYg)  =  i 


-30(4+50+602+03)/[(l+0)2(2+0)2] 
-3(1-0) (2+30-02)/[(l+0)2(2+0)] 


°  <  6  <  h  , 
H  <  B  <  1  . 


(4.9) 


A  graph  of  this  correlation  function  is  given  in  Figure  2;  its  maximum 
negative  value  is  -13/15  =  -0.8667  which  occurs  at  p  =  0.5;  this  indicates  for 
the  EP2  model  a  good  negative  range  to  the  Spearman  correlation,  but  no 
ultimate  degeneracy.  Bivariate  scatters  from  the  UP2  model  show  a  strong 
degree  of  linearity  compared  to  the  near  antithetically  curved  scatter  of 
the  corresponding  exponential  EP(2)  pair  in  Figure  5.  Again,  the  asymmetry  of 
the  dependency  about  6  =  0.5  in  the  EP2  model  is  more  evident  from  the 
Spearman  correlation  of  the  UP2  model.  Negative  dependency  is  considerably 
stronger  for  6  *  0.5  -  <5  than  the  corresponding  0  =  0.5  +  6  .  This  parallels 
the  EP1  model  in  which  positive  dependency  is  stronger  for  8  =  0.5  -  6  than 
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for  3  =  0.5  +  6  ;  these  remarks  are  borne  out  in  scatters  for  various  values 
of  3  •  Finally,  it  should  be  noted  that  the  antithetic  choice  (I,a(I))  is 
best  in  regard  to  obtaining  negative  dependency  in  the  UP2  model,  as  was  found 
for  the  EP2  model.  Regressions  are  again  not  easily  tractable. 

Positive  and  Negative  Dependence 

Extension  of  the  bivariate  Models  given  above  to  models  which  exhibit 
both  positive  and  negative  product-moment  correlation  is  achieved  quite 
simply  by  mixture  arguments.  For  example,  combining  (2.3)  and  (4.2)  let 

Xj  =  BEj  +  I-jEg 

f BE2  +  1}  E-j  w.p.  p  ,  (4.10) 


(3a(E2)  +  I2a(E1)  w.p.  1-p  . 

This  uses  the  fact  that  a  probability  mixture  of  two  identically  distributed 
random  variables  (possibly  dependent)  has  that  same  identical  distribution. 
These  two  parameter  model  contains  both  the  EP1  model  (p  =  1)  and  the  EP2 
model  (p  =  0)  and  has  product  moment  correlation  ranging  between  -0.5724 
and  0.75. 
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5.  TWO  PARAMETER  VERSIONS  OF  THE  MODELS:  EP3,  UP3,  EP4,  UP4  MODELS 
Dependency  in  the  models  has  thus  far  been  determined  by  one  parameter 
6  ,  but  it  may  be  desirable  to  have  an  additional  parameter  in  the  model  so 
that  two  aspects  of  the  joint  distribution  can  be  modelled;  for  instance, 
one  dependency  measure,  such  as  the  product-moment  correlation,  and  a  prob¬ 
ability  statement,  such  as  P(X1  >  X2),  which  in  the  EP1  model  was  equal  to 
0.5  for  all  g  .  Two  approaches  will  be  discussed  rather  briefly.  One 
allows  6  in  the  randomly  linear  combination  of  E-|  and  E2  to  be  replaced 
by  g-|  for  X-j  and  g2  for  X2  ,  while  in  Section  6  a  two  parameter  randomly- 
1  inear  operation  is  used:  this  stems  from  the  NEAR(l)  model  of  Lawrance 
[1980],  Lawrance  and  Lewis  [1981]. 

The  EP3  Model 

The  first  two  parameter  model,  EP3,  is  defined  by  the  equations 

X-j  =  g-j  E-j  +  I  ■j  E2  ,  0  g-j  <_  1  , 

(EP3)  (5.1) 

x2  =  g2  e2  1 2  E i  ,  0  b2  5.  ^  • 

The  marginal  distributions  of  1^  and  I2  must  correspond  to  the  different 
g  parameters  so  that  P(I-j  =  0}  =  1  -  P{  1^  =  1}  =  g^  and 
P(I2  =  0}  =  1  -  P{I2  =  1}  =  g2  .  The  product  moment  correlation  of  the 
model  is  given  by 

Corr(XrX2)  =  B1(l-B2)  +  g^l-g^  +  Covd-,  ,I2)  .  (5.2) 

The  general  joint  distribution  for  (I-|,I2)  is  the  obvious  generalization  of 
(4.4);  the  value  of  p^  for  maximum  positive  dependence  is  now  min(l-g^ ,l-g2) 
Hence  the  result 


Corr(X1,X2)  =  < 


[26,(1-62)  +  62(l-6,) 

6,0-62)  +  2620-6,) 


0  <.  3,  <  62  1  1  » 
0  £  62  £  B,  £  1  • 


Contours  of  this  function  (Figure  6}  show  a  central  region  with  a  local 
maximum  of  0.75;  the  function  decreasing  from  each  side  0*  the  centre  along 
the  line  6,  *  62  to  zero  values  at  the  corners  where  6,  =  =  0  and 

61  =  62  =  1  »  and  the  function  increasing  from  each  side  of  the  centre  along 

the  line  6,  +  62  =  *  to  unit  values  at  the  other  two  corners.  The  joint 

density  of  X, ,X2  is  simply  obtainable  but  is  now  a  mixture  of  three, 

instead  of  two,  components.  Thus  maximum  likelihood  estimates  of  6,  and  e2 
can  be  obtained. 


Figure  6.  Correlation  function  of  the  two-parameter  positively 
dependent  EP3(6, ,02)  model:  unlabelled  contours  may  be  identified 
from  their  join!  wfth  the  6,  and  62  axes. 
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The  edges  (6^  =  0)  or  (B£  “  0)  give  a  one-parameter  bivariate  exponen¬ 
tial  pair  which  corresponds  to  adjacent  (exponential)  variables  in  an 
EAR(l)  process  [Gaver  and  Lewis,  1980].  Note  that  using  the  basic  relation¬ 
ship  (2.1)  outside  of  the  serial  constraint  posed  by  a  time  series  gives  a 
broader  bivariate  exponential  pair. 


The  UP3  Model 

The  corresponding  uniform  model,  UP3,  is  given  by 


(UP3) 


^1  ^ ?  ^2 

Y,  -  1  -  V^1,  Y2  =  1  - 


(5.4) 


and  has  correlation  (the  Spearman  correlation  of  the  EP3  model),  given  by 


Corr(Y1,Y2)  =  12 


1-6, 


e2'el 


(2+e1 )(i+e2)  +  (2+82j(i+81) 


(5.5) 


+  (H^HT-^) 


-  3  (o  <  ^  <  a2  <  i ) 


For  6^  1  ^2  there  1S  a  symmetrical  result.  The  contours  of  this  function 
represent  a  deformation  of  those  just  described  for  the  EP3  model,  symmetry 
about  6}  =  82  is  preserved,  while  the  area  of  local  maximum  is  moved 
towards  the  origin.  Thus,  asymmetry  in  the  Spearman  dependence  is  indicated. 


The  EP4  Model 

The  negatively  dependent  version  of  the  EP3  model  is  the  two  parameter 
version  of  (4.2),  and  will  be  called  the  EP4  model;  this  has  correlation 
function 
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I 


i 


Corr(X],X2)  =  -  0.64491^ (l-02)  +  ^(1-^)] 

@1^2  for  t  62  £  1  , 

—  • 

(l-31)(l-$2)  for  +  32  1  1  • 


(5.6) 


Contours  of  this  function  indicate  a  valley  of  lowest  negative  values  along 
the  direction  8-]  +  82  =  1  ,  with  a  central  value  of  -0.5720  and  corner  values 
of  maximum  negativity,  -0.6449.  For  specified  negative  correlation  down  to 
about  -0.5  there  is  a  wide  choice  of  (B-|  »B2)  pairs,  in  the  two  regions  either 
side  of  the  valley  along  8-j  +  82  s  1  • 


The  UP4  Model 

The  final  development  of  the  two  parameter  models  is  to  the  negatively 
dependent  two  parameter  uniform  model,  the  UP4  model.  This  corresponds  to 
(4.8)  with  8  and  I  in  Y1  replaced  by  8-j  and  I-j  and  8  and  I  in 
replaced  by  82  and  a(I2)  .  Maximum  negativity  of  dependence  follows 
from  the  'antithetic'  joint  distribution  for  1^  and  I2  which  has 
P-jj  =1  -  (8^+B2)  when  8^  +  82  1  1  and  pQ0  =  8^  +  82  -  1  when 
+  $2  —  ^  '  We  ^en  ^ave»  f°r  t^ie  case  6-j  +  82  £  1  » 


Corr(Ui ,U2)  =  12 


1  -  (3-,  +62) 


(l+81)(2+81T(l+82)(2+82)  + 


+  (l+81)(l+82)(2+82) 


(5.7) 


The  corresponding  result  for  8-j  +  82  £  1  is  omitted.  The  contours  of  this 
function  are  quite  similar  to  those  of  the  EP4  model,  but  with  a  central 
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value  of  -0.867  and  associated  corner  values  of  minus  one;  the  other  two 
corner  values  are  zero,  as  for  the  EP4  model;  it  is  again  useful  to  view 
these  as  the  Spearman  correlations  of  the  EP4  model. 

From  the  point  of  view  of  simulation,  the  two  parameter  models  are 
hardly  more  complicated  than  the  one  parameter  version  and  offer  an  extra 
degree  of  flexibility.  The  joint  probability  density  functions  and  moment 
generating  functions  are  more  complicated  than  for  the  one  parameter  models, 
but  can  be  derived  if  needed,  for  instance,  for  estimation.  It  is  also 
possible  to  derive  measures  such  as  P{X^  >  Xg)  ,  which  is  not  necessarily 
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6.  MORE  GENERAL  FOUR  PARAMETER  MODELS:  EP+,  EP-MODELS 


The  NEAR(l)  exponential  time  series  models  of  Lawrance  and  Lewis  [1981] 
suggests  a  class  of  four-parameter  models;  these  are  likely  to  be  over¬ 
complicated  for  ordinary  use.  but  place  the  earlier  models  in  a  general 
setting  and  suggest  a  further  class  of  two-parameter  models.  They  also  give 
a  bivariate  exponential  pair  with  full  positive  correlation,  and  without  the 
degeneracy  which  occurs  for  the  case  (8^  *  0  or  e2  *  0)  in  the  EP3  model. 
Being  as  brief  as  possible,  this  development  gives  the  EP+  model  as 


(EP+) 

where,  for  i  =  1,2, 

1  w.p.  cu 
Vi 

0  w.p.  l-ai 


X1  =  31V1E1  +  !1E2  ’ 

(6.1) 

X2  =  ^2^2E2  +  *2E1  * 

1  w.p.  (l-f^/Cl-O-a.)^] 

*i  M 

(1-a..  )|3.j  w.p.  ai81./[l-(l-ai)Bi] 


and  the  random  variables  (V-j , V2) ,  (I-|,I2)  are  independent  between  pairs  but 
usually  dependent  within  pairs.  This  model  has  correlation  structure  of 
the  form 


Corr(XltX2)  =  8^  CovO^.Vg)  +  Cov(I1 ,12) 

(6.2) 

^  ®1 6)  ( 1  “^2^2^  a2^2^"otl^l  )  * 

For  a  model  of  maximum  positive  dependency  the  appropriate  joint  distributions 
of  (V-|,V2)  and  (Ij,I2)  could  be  obtained  using  constructions  as  at  (4.4). 

A  further  two  parameter  model,  EP5,  is  suggested  by  taking  ot^  =  a2  =  a  , 
81  «  02  »  8  and  V1  =  V2,  ^  *  I2  ;  then  (6.2)  reduces  to 
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Co rrU-j.Xg)  =  3ctB(1-a6)  . 


(6.3) 


This  has  maximum  value  0.75  at  ctB  *  0.5  ,  but  does  not  attain  the  higher 
correlations  of  the  EP3  model,  and  Is  more  complicated  In  Its  construction. 

Both  the  EP1  and  EP3  models  can  be  obtained  as  special  cases  of  the 
EP+  model;  the  EP1  model  is  given  by  taking  =  ot2  =  1  and  B-j  =  &2  =  ^ 
with  1^  =  I2  .  while  the  EP3  model  has  a-j  =  a2  =  1 ,  (6-j *B2)  unchanged,  and 
with  (I-j,I2)  of  maximum  possible  dependency. 

Corresponding  to  (6.1)  there  is  a  more  general  negatively  dependent 
four-parameter  model,  EP-,  of  the  form 


X1  =  61V1E1  +  T1E2  * 

X2  *  B1V1a(E2)  +  , 


(6.4) 


which  has,  by  analogy,  the  correlation  structure 


Corr(X-|  ,X2)  *  B1B2Cov(V1  ,V2)  +  Cov(IrI2) 

-  0.6449[a1B-|(l-a2^2^  +  a2e2^-al^l^  ' 


(6.5) 


It  contains  the  EP2  and  EP4  models,  and  would  suggest  the  further  two  param¬ 
eter  model,  EP6,  by  taking  =  a2  ■  a,  Bj  *  B2  *  6  and  (V^,V2),  (Ij,I2) 
as  two  antithetic  pairs. 

Any  further  details  are  omitted,  as  are  the  corresponding  uniform  models 
of  the  four- parameter  class,  UP+,  UP-  given  in  Table  1,  which  summarizes  all 
the  models. 
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7.  MULTIVARIATE  GENERALIZATIONS 

The  possibility  of  multivariate  generalizations  is  apparent  by  repeated 
use  of  the  random-coefficient  linear  functions  of  (2.1).  In  fact  this  is 
the  way  the  EAR(l)  process  is  constructea  with  a  serial  chaining  of  i.i.d. 
exponentials  E-j ,  E2,....  Any  k  variables  in  an  EAR(l)  process  (or  an 
NEAR(l)  process)  are  a  k-variate  exponential  random  variable.  Outside  of 
this  serial  context  many  possibilities  suggest  themselves.  Thus  changing 
9  to  gj  ,  I  to  1}  and  replacing  E2  by  a  similar  random-coefficient 
linear  function  of  independent  E2  and  E^  using  e2  and  I2  ,  gives 

X  =  E-j  +  6 2 T -j E 2  ^  *  (7.1) 

Here  1^  and  I2  must  be  independent.  A  triple  of  dependent  exponential 
variables  could  be  constructed  using  E-j ,  E2,  E3  in  three  of  the  six 
possible  orders. 

A  simpler  possibility  would  be  to  use  the  basic  randomly- linear 
operation  on  the  pairs  (EpE2)  ,  (E2>E3)  and  (E^.E^)  .  Such  developments 
are  not  considered  here.  The  possibilities  are  legion  but  no  'natural' 
simple  method  suggests  itself  above  any  other. 
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1 

I ' 

I  1 

8.  LINE  DISCONTINUITIES  AND  RELIABILITY  PROBLEMS 

r 

The  problem  of  discriminating  amongst  the  many  bivariate  exponential  , 

models  is  not  simple,  although  there  may  be  modelling  or  structural  details 

* 

which  recommend  certain  models  in  certain  contexts.  Most  of  these  details  } 

are,  however,  tenuous  and  difficult  to  verify  from  data.  An  important 

mathematical  scheme  is  given  in  Griffiths  [1969],  who  found  a  canonical  ! 

expansion  for  bivariate  exponential  random  variables.  Again  the  models  put  | 

forward  here  are  simple  to  generate  on  a  computer  and  are  analytically  ’! 

tractable.  ; 

Another  property  is  that  the  bivariate  distributions  do  not  have  line 

1  4 

discontinuities.  However,  as  an  illustration  of  modelling  considerations, 

i 

we  note  that  in  a  reliability  context  thus  is  not  necessarily  a  virtue.  It 
is  known  that  components  in  a  system  can  fail  from  a  common  cause,  which  is 
precisely  what  gives  the  line  discontinuity  in  the  bivariate  distribution. 

* 

However,  it  is  simple  to  put  this  in  to  the  present  models  and  it  can  be 

done  in  at  least  three  ways.  j 

Thus,  let  (XpX2)  denote  any  unit  mean,  bivariate  exponential  pair, 
let  (Ipl2)  denote  an  indicator  pair,  possibly  completely  or  partially  , 

dependent,  with  marginal  distributions  P{I^  *  1}  =  1  -  P{I^  =  0}  =  P1  and 

;  i 

P{ I2  =  1}  =  1  -  P{l£  =  0}  =  >  ar,d  E  be  an  independent,  unit  mean 

exponential  random  variable.  Three  new  bivariate  exponential  pairs  are  given  , 


(1) 

Z1  =  I]X1  +  (1-I-j )E  , 

Z2  =  I2X2  +  (1-I2)E  ; 

(8.1) 

(11) 

Z1  *  BE  +  I1X]  , 

Z2  =  BE  +  I£X2  ; 

(8.2) 

and 
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( 1  i  1 )  Z-j  =  min(X.j,E)  ,  Z2  =  min(X2,E)  l  .3) 

In  all  three  cases  there  is  a  non-zero  probability  that  Z^  and  Z2  are 
proportional  to  E  . 

The  first  pair  uses  the  idea  that  mixtures  of  identically  distributed 
random  variables  have  that  same  distribution,  the  second  pair  uses  the  basi 
relationship  (2.1)  and  the  third  uses  the  fact  that  the  minimum  of  Independent 
exponential  random  variables  is  an  exponential  random  variable.  In  fact, 
if  Xj  and  X2  are  independent,  then  (8.3)  is  the  Marshal  1-01  kin  model. 
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TABLE  I 

Summary  of  Models  Considered  {See  text  for  relevant  details) 


EP1 

X1  =  BE1  +  IE2 

X2  =  BE2  +  IE1 

UP1 

Y1  =  1-U^U| 

y2  =  i-u{ujj 

EP2 

X1  =  *1  +  ¥2 

X2  =  3a (Eg)  +  1 ( E -j ) 

UP2 

Y]  =  1-UjuJ 

y2  =  l-(l-U1)a<I>(l-U2)e 

EP3 

X1  =  B1E1  +  ¥2 

X2  =  ^2E2  +  ¥l 

UP3 

S-I  r 

V,  - 

1  p  Bp 

Y  =  l-U  d 

2  1  1  2 

EP4 

X1  =  ¥l  +  ¥2 

X2  *  ( ^2 ^  +  ¥^E1 ^ 

UP4 

I 

Y1  =  ^l  U2 

Y2  -  l-O-u/^i-u,/2 

EP+ 

X1  =  61V1E1  +  ¥2 

X2  =  &1V2E2  +  ¥l 

UP+ 

^1^1  *1 

Y-j  =  1-U^  'U2' 

IP  BpVp 
y2  =  i-u^u2  c 

EP- 

X1  =  31V1E1  +  ¥2 

X2  =  B^afEg)  +  I2a(E1) 

UP- 

^i v i  I, 

^  *  l-u/  '(1-U2)  1 

Ip  3-iVp 

y2  *  l-O-u,)  2(1-U2)  1  2 
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